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Abstract. 

We analyze the role of impurities in the fractional quantum Hall effect using a highly 
controllable system of ultracold atoms. We investigate the mechanism responsible for the 
formation of plateaux in the resistivity/conductivity as a function of the applied magnetic 
field in the lowest Landau level regime. To this aim, we consider an impurity immersed in 
a small cloud of an ultracold quantum Bose gas subjected to an artificial magnetic field. We 
consider scenarios corresponding to experimentally realistic systems with gauge fields induced 
by rotation of the trapping parabolic potential. Systems of this kind are adequate to simulate 
quantum Hall effects in ultracold atom setups. We use exact diagonalization for few atoms and, 
to emulate transport equations, we analyze the time evolution of the system under a periodic 
perturbation. We provide a theoretical proposal to detect the up-to-now elusive presence of 
strongly correlated states related to fractional filling factors in the context of ultracold atoms. 
We analyze the conditions under which these strongly correlated states are associated with the 
presence of the resistivity/conductivity plateaux. Our main result is the presence of a plateau 
in a region, where the transfer between localized and non-localized particles takes place, as a 
necessary condition to maintain a constant value of the resistivity/conductivity as the magnetic 
field increases. 


PACS numbers: 03.75.Hh, 03.75.Kk, 67.40. Vs 


Quantum simulation of conductivity plateaux and fractional quantum Hall effect ... 2 

1. Introduction 

Transport properties provide some of the most fundamental characteristics of condensed 
matter systems (cf. []Tl|2]|). In contrast, in physics of ultracold atomic and molecular gases 
0, the studies of transport, unlike those performed in solid state settings, are hindered by the 
difficulty of having continuous and durable flow of atoms; for this reason they have been very 
limited so far. Among others they included: the investigations of Bloch oscillations (from 
the early studies with cold atoms 0 to the recent experiments with disordered gases 0), the 
extensive work on transport and diffusion in disordered gases [[6l EL S IH [[Oil , and the very 
recent experiments on quantized conductivity ifTTT fl2l [T31IT41 15]. 

Paradigmatic systems, in which the transport properties play an essential role, are the 
systems that exhibit integer or fractional quantum Hall effects (IQHE, or FQHE) [fl6l 117]. 
The quantum Hall effect consists in fact in quantisation of the transverse conductance for 
electronic current in the condensed matter systems, and for atomic flow for neutral atomic 
gases. Although in the IQHE the interactions play an irrelevant role, the underlying physics, 
even if well understood, is highly non-trivial. The case of the FQHE, where the interaction has 
a crucial contribution, is more complex and not yet completely understood. For these reasons 
systems that exhibit FQHE belong to the most popular systems of strongly correlated particles 
that still await conclusive explanations and “call for” quantum simulations, for instance with 
ultracold atoms or ions 0. 

In order to quantum simulate QHE it is necessary to generate strong artificial magnetic 
(gauge) fields. In the context of ultracold atomic systems, first realizations of artificial 
gauge (magnetic) fields were considered in rotating traps [HIQ3L 20;[2H 22, 23, l24ll . Quite 
soon it was realized, however, that the most promising way to generate the artificial gauge 
fields is to use the laser induced fields - these methods are described in detail in several 
reviews 0i25 ,;26 ; ;|27], while for the recent spectacular experiments the reader should consult 
cf. [28 , [29, 30, 3T1 . Various methods of detection of the Hall effect have been proposed and 
realized in the mentioned experiments (see cf. [321 1331). In particular, it was shown lf34l 
how the quantized Hall conductance can be measured from density profiles using the Streda 
formula lf35fl . In this paper we propose to measure the quantized Hall conductance directly 
as a transport property and suggest to use the response of the considered system to the time 
dependent perturbation. 

Before we turn to atomic systems, it is instructive to review briefly the phenomenology 
of electronic systems [16l [36]]. Hall effect appears already in classical physics, where the 
transverse resistivity is proportional to the magnetic field B. The transverse conductivity may 
be expressed by the famous expression 

e 2 

CTyx = V—, ( 1 ) 

where v = nh/eB is the filling factor, and n is the electron areal density. Amazingly the 
same formula holds in the quantum mechanical case, being the consequence of the Galilean 
invariance ll36ll . 

The explanation of plateaux corresponding to integer filling factors observed in 
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experiments requires thus additional arguments. These arguments are based on the fact that 
in typical experimental situations Galilean invariance breaks down due to the presence of 
random impurities. Accordingly, the spectrum of (non-interacting) 2D electron gas in the 
magnetic field does not exhibit discrete Landau levels only. Close to the Landau levels in 
fact the spectrum consists of a band corresponding to extended (conducting) states. Far from 
the Landau level energy the spectrum corresponds to states localized due to the presence 
of impurities via the mechanism of Anderson localization 0711 . Obviously, localized states 
do not contribute to the conductivity, and thus one can expect that when Fermi energy 
decreases between two subsequent Landau levels including less and less localized states 
in the Fermi sea, keeping the condition of fully occupied Landau levels, the conductivity 

does not change i.e. it exhibits a plateau. Why has the plateau the value exactly equal to 

2 

c>yx = v%- with v integer, and why is this result so robust was a rather surprising fact in 
the beginning of the 1980s. It was first explained by the famous Laughlin argument lf38ll . 
He demonstrated the quantization of the Hall resistivity analyzing an imaginary experiment 
topologically equivalent to a ”Corbino-type” sample of a disk shape with a central hole llT6ll . 
Laughlin’s arguments were then generalized by Halperin |f39l and Biittiker ll40ll to the strip 
geometry, employing the properties of the edge states and edge currents. Contemporary 
understanding of the robustness of the IQHE is based on the topological nature of integer 
transverse conductivity, first related to Chern numbers by Niu et al. [14TI (see also Ii42l and 
references therein). 

The IQHE requires high, but not extensively high values of the magnetic field, in 
which several Landau levels are involved. The quantization of the transverse conductivity 
corresponding to integer values of the filling factor, is of course due to the quantization of the 
Landau levels. Still, the step-wise behavior of conductivity in the 2D electron gas originates 
from the influence of impurities ll43l . 

In contrast, in the case of strong magnetic fields in the lowest Landau level regime, the 
non-interacting particle approach cannot be applied. One can think of the composite-fermion 
picture in which the fractional filling factor for electrons is transformed into integer filling 
factors for new quasiparticles: electrons dressed with magnetic quantum fluxes fill completely 
several Landau levels [44, 45J. However, for this equivalent system of composite fermions in 
the IQH regime, the role of impurities is strongly combined with effects of interaction. One 
can think that the impurities play the role of a reservoir of particles trapping or releasing 
particles as the Lermi level moves across the localized states as the magnetic field changes. 
Notice that the Lermi energy decreases to lower values as the real magnetic field B grows 
m. As a consequence, for some intervals of B the density of the extended electrons, those 
that contribute to the current, increases due to the transfer of electrons from impurities to 
the Landau levels. This effect compensates the increase of B, providing the appearance of a 
plateau in the resistivity/conductivity: 

p yx ~ B/n e , (2) 

where n e is the density of the extended part of the system. It must be stressed that the 
presence of impurities plays the same role as that of the edge in finite systems. The main 
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ingredient is the presence of a scalar potential locally linear in x that traps the particles Ifl6l . 

The size of the plateaux depends thus on the number and the properties of the impurities. 
It must be realized that these plateaux appear on special values of p yx that localize states 
of significant interest, with fractional filling factors. Without impurities these values of the 
resistivity would not be visible. 

Turning back to the atomic gases, to simulate the similar phenomenology, we must 
somehow obtain localized and non-localized particles and look for regimes where transfer 
between them is possible. Since in ultracold atom setups one has the possibility to engineer 
controllable impurities, such systems provide an ideal tool to understand the role of impurities 
in the formation of the plateaux and their interplay with interactions in the FQHE. In our 
numerical simulations, the possibility of distinguishing localized and non-localized particles 
was achieved in the following way: The diagonalization of the one-body density matrix 
provides us with the natural orbitals. Importantly, in all the analyzed cases, one of the orbitals 
is mostly concentrated around the impurity. In contrast, the other ones remain extended. 
Therefore we can distinguish between these two parts. Intervals of the artificial magnetic field 
B* where the occupations of the natural orbits have a significant variation with B* turned out 
to be crucial to identify the regions where transfer is possible and plateaux are expected. 

Recently we have used state-of-art exact diagonalization to study properties of small 
clouds of atoms in a trap under influence of strong artificial gauge fields (see |[22l 23, 46, 47. 
148 , 49, SOlEjj). In this paper, we expand the previous studies [|52l . Our main goal is to learn 
about the relationship between the presence of impurities and the appearance of plateaux in 
the Hall resistivity as a function of the magnetic field in the fractional quantum Hall regime. 
This regime has not been achieved experimentally with atoms and our preliminary analysis 
is intended to predict possible future results. We use here the exact diagonalization method 
to calculate the ground state (GS) and its excitations in the absence/presence of an impurity. 
We analyze the time evolution (TE) of the system submitted to a periodic perturbation which 
represents an applied external electric field to simulate the transport equation 

jy (^) t7 yx E x (lfj ; (3) 

where j y is the equivalent to the electronic current for atoms and o yx is the transverse 
conductivity. In the Appendix we compare the results with those obtained using the linear 
response theory (LRT) approximation and conclude that, aside some limitations associated 
with resonances, the comparison is extremely good. 

Our main result is the appearance of a plateau close to the GS transition that takes place 
when the angular momentum changes from L = N(N — 2) to L = iV(iV — 1). The change 
of angular momentum modifies the resistivity and produces a bump partially overlapped with 
the plateau. An important outcome is that the presence of an impurity is a necessary condition 
to generate plateaux. With no impurities the change of the occupations is abrupt and the 
transfer process is not possible. 

The numerical complexity of the problem allows us to study only rather small systems 
up to N = 4 atoms. In effect the predicted plateau is small. We expect, however, that for 
large systems (N = 100-1000) the natural increase of the number of impurities and thus the 
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increase of the localized part, would guarantee the visibility of the plateau. On the other hand, 
the robustness comes from the topological nature of the conductivity, as previously mentioned. 

This paper is organized as follows: In section II we present the model for the basic 
Hamiltonian H 0 to calculate the full spectrum of the unperturbed system. From its GS, we 
obtain the natural orbitals and their occupations. In addition, we show the explicit expression 
of the periodic perturbation and the expression for the equivalent periodic trapping potential in 
the full Hamiltonian. In section III, we analyze the time evolution of the expected value of the 
current operator j y and identify the conductivity. The main results are shown. In section IV 
we define a kind of restricted operators that emphasize the presence of the plateau and discuss 
about its meaning. Finally section V contains the summary and discussion. The appendix 
develops the LRT and shows some comparisons with the TE, as well as some limitations 
inherently related with the method. 


2. The model 


We consider a system of N one-component bosonic atoms of mass M confined on the XY- 
plane. The cloud is trapped by a rotating parabolic potential of frequency u±_ and rotation 
Q along the Z-axis, rotation which in an effective way generates an artificial magnetic field 
denoted by B*. In the rotating reference frame the basic Hamiltonian (not including the 
perturbation) reads 


Hq H sp T Hint: 

where the single particle part is given by 

TT \" r P M 2-2 
H sp = 2^ \ 7TTT + — W | r 


i= 1 


12M 


Hz - L 


+ W 


which can be rewritten in an equivalent way as 


N 

4p= £ 

i =1 


— (p + A) 2 + -M (ul- ll "‘ 
12M yv ’ 2 I x 


* \ 2 


AM 2 


+ W, 


with 


B* - B* 

A = -ttV , Ay = -—X. 


(4) 


(5) 


( 6 ) 


(7) 


The particular selection of the symmetric gauge has been done in the definition of A, with 

B* = 2 Mfi (8) 

being the modulus of a constant artificial magnetic field directed downward along the Z- 
direction and r = (x, y). 

The potential W = ff 1 . , Wj is due to the presence of K impurities, which are modeled 
by Dirac delta functions, 


h 2 


N 


Wj = 


i= 1 


(9) 
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Figure 1. (Color online) Density of the first five natural orbitals, the only ones with non- 
negligible occupations for N = 3 (upper row) and for N = 4 (lower row). For N = 3, 
the second one, panel (b), is localized close to the impurity siting at ai = (1,0). In 
contrast, for N = 4, the natural orbital which is localized around the impurity is the first 
one [panel (f)]. The occupations are: 0.92,0.74,0.57,0.42 and 0.34 for N = 3 and 
1.30,0.83,0.68,0.52,0.41 and 0.26 for N = 4, respectively. In both cases we take Ng = 6, 
71 = 0.1, Q = 1.95. 

The dimensionless parameters 7 y measure the strength of the impurities and aj localize them 
on the XY plane. The term W breaks the circular symmetry except for the case of a single 
impurity localized exactly at the center. 

We model the atomic interaction by a 2D contact potential characterized by 

Hint » ( 10 ) 

i<j 

where g = \/&na/\ z is the dimensionless coupling, a is the 3D scattering length and 
A, = \Jk/Muj z . We assume oj z , the trap frequency in the Z-direction, much larger than any 
of the energy scales involved, in such a way that only the lowest level is occupied. Therefore, 
the dynamic is frozen in the Z-axis and the system can be considered two-dimensional. 

Without impurities, the solutions of the single particle part produce the Landau level 
structure If45l . The energy levels are separated by h(u>± + Q). We assume that B* is large 
enough to consider just the lowest Landau level (LLL) regime where the appearance of energy 
gaps has a completely different origin as those in the IQH, where several Landau levels are 
implied. Within this regime, the kinetic part of the Hamiltonian reads 

Hkin = h(u± — Q) L + Nhw± . (11) 

The single particle solutions with well defined angular momentum m are the Fock-Darwin 
(FD) functions, given by (f> m (6,r) = e~ r2 / 2 r m . 

Once the spectrum of the whole system is obtained for a given number of particles N 
and for fixed values of Q, aj and g, we proceed to distinguish the eigenfunctions of 
the one-body density matrix according to their localization. To this end, we diagonalize the 
one-body density matrix given by 

p ( 1 ) (r,r') = (# t (r)^(r / )), 


( 12 ) 
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where the expected value is calculated at the GS and 'h(r) is the field operator. Its 
eigenfunctions are the natural orbitals ifi , linear combinations of the FD functions, and the 
eigenvalues are their occupations n*, i = 1 , l m + 1 (l m varied until convergence). 

From now on, we use the complete set of natural orbitals as a base to represent functions 
and operators in the second quantized formalism. Next we analyze the density distribution 
of each orbital and look for their localization around the impurities. For N = 3 and one 
impurity at a! = (1,0) [see Eq. ([9]) - lengths are in units of the XY harmonic oscillator 
length], the result is that the orbital out of five orbitals with non-negligible occupation 
(the orbitals are ordered by decreasing occupations), presents a density distribution mainly 
localized at the impurity, as is shown in Fig. |T| (upper row). Similar results are obtained for 
N = 4 where the most localized orbital is Vh (lower row). This localization of some orbitals 
allows us to distinguish between localized and extended states. 

Finally, we consider that the cloud of atoms is dynamically forced by an oscillating term, 
while the impurity remains attached to a fixed position. To be explicit, the full Hamiltonian 
of the system is 

H{t) = Hq + H P ert(t) (13) 


Hpert(t) = -A f y^ y Xj J £{t) sin (ut) = ^2 f(t) (I 4 ) 

\ i=l J *=1 

where A gives the intensity of the perturbation, which we assume small. The explicit form of 

£(t) is 

£{t) = 1 - exp[-(t/cr) 2 ] (15) 

where a determines the velocity of the evolution. The perturbation is switched on at t = 0 
and, as t increases, the stationary regime is achieved when the amplitude of the oscillations 
can be considered constant. From now on we consider M — 1/2 and h = 1 and choose 
Aj_ = \Jjt~ = \/2M, huJ±/2 and oo±/2 as units of length, energy, and frequency, 
respectively. With our unit of length, lu± = 2. In the simulation, the sequence — A£(£) sin (ut) 
is identified with the electric field E x (t) in the transport equation. Namely, for a single particle 
and a single impurity [see Eq. ©], including the perturbation we have 

Hi = (p + A)- + {u]_ - f) 2 )(x 2 + y 2 )i + f{t)±i - 7i5 (2) (fi - ai). (16) 


The effective trapping potential can be re-written as an oscillating trap 

.2 r>2\ T is-. , /(^) 


(wI-H 2 ) (x + 


2(cu 2 - H 2 ) 


) 2 + y 2 


(17) 


3. Time evolution 

Let us show that this model allows us to identify the transverse conductivity from the transport 
equation £2). For that, we need to analyze the time evolution of the expected value of the 
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Figure 2. (a) Occupations of the first two orbitals as a function of O. There is a correlation 
between the steps of the occupations and the jumps of (L) in a nearly symmetric system. For 
N = 4 the magic numbers of L in the symmetric case are 0 — 4 — 8—12. The vertical lines 
at 1.765 ,1.920 and 1.959 mark the position of the jumps of L for the symmetric case. We 
considered N = 4, Ng = 6 , oj = 1.1 71 = 0.1 and ai = (1,0). (b) A zoom in the 
interesting region [which is highlighted with a small rectangle in (a)]. The localized orbital for 
N = 4 is fi. 


current operator j y , which is given by 

1 


jy( r ) = [Py + A y (r)J ¥(r) 

Py~M r )] ^( r )l- 


(18) 


By calculating (j y ) t once the stationary regime is reached and in the case that we obtain a 
linear behavior in A [see Eq. ([T4l) l. we are able to obtain the transverse conductivity from the 
transport equation © due to the identification we have done between the perturbation and the 
electric potential associated with a constant electric field. 

To obtain (t)\j y \^(t)) we solve the Schrodinger equation id t ^{t) = H(t)^(t) 
with the time-dependent Hamiltonian given by Eq. (fl3l) . We consider the wave function as 
T(f) = J2n =1 c n(t) where the set {$„.}, n = 1, is a basis of the Hilbert space 
of dimension rid, given by the many-body wave functions which solve Hamiltonian H 0 with 
eigenvalues {E n }. Then, we obtain the system of equations: 


tdt Cn(f) Cn(t)E n 

— A f{t)sin(ujf) Ec m (t) < n\x\m > 


nd 


(19) 


m= 1 


with the GS($i) as the initial condition, i.e., ^(f = 0) = $ 1 . We solve these equations using 
the Runge-Kutta fourth-order algorithm. Once we obtain the transverse conductivity o yx , we 
can obtain the resistivity from 


cr 


Pyx 


yx 


<T 


yx 1 


+ Wx 


( 20 ) 


where o xx is obtained from j x = a xx E x . 

Let us note here that our evolution is rather adiabatic. Indeed, we have checked that the 
overlap |('T G s|T'(t))| remains nearly one at all times. E.g., we numerically obtain that, for 
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Figure 3. (a) Hall resistivity p yx as a function of Q. We used N = 4, Ng = 6 , uj = 1.1, 
71 = 0.1, A = 0.001. We considered an impurity at the X-axis at different positions, that 
is, ai = (xo, 0). The vertical dashed lines mark the interval considered in panel (b), where 
we show a zoom of the area of interest. The vertical dotted line at 1.962 in panel (b) marks 
the approximate place of the jump in the angular momentum from 8 to 12. We also include 
two additional positions of the impurity to show that the effect depends on the position of the 
impurity. 


o = 10 (see Eq. (fl5l)j. 

|(tf GS |tf(tc)>| 2 = 0.97 (21) 

where 'kgs is the ground state of H(t c ) — H 0 + H pert (t c ) and \k(£ c ) is the solution of the 
Schrodinger equation at t c (chosen such that sin[ut c ) = 1). 

In agreement with our previous discussion, to generate a plateau we must look for an 
interval of where the occupation of the localized orbital changes. To this end, we analyze 
the orbital occupations as a function of [see Fig.[2](a)] and focus on the region fulfilling two 
requirements: On one hand, the occupation of the localized orbital decreases as B* increases 
producing an increase of the extended part and, on the other hand, this region lies within the 
largest possible value of (L) where plateaux are expected. In other words, there are two kind 
of intervals: i) intervals where the localized density is flat giving linear dependence of p yx 
with B* where quantum and classical behavior coincide lfl6l : and ii) intervals where a plateau 
occurs and then, the change of B* drags p d (the extended density), and p yx remains constant. 
Notice that this means that for our analyzed small samples we expect only one plateau along 
the whole interval of the largest value of (L ). Fig. |2](b) is a zoom of the interesting region. 

Our main result is represented in Fig. [3j where we show the variation of the Hall 
resistivity as a function of Q for some values of the position of the impurity, which we 
always assume to be in the A"-axis, that is a! = (x 0 , 0). In Fig. [3] (a) we show the Hall 
resistivity for three exemplary values of the position of the impurity for the whole range of H 
investigated. We observe that some structure, different from the linear behavior, occurs just 
in the region shown in Fig.[2](b). This structure begins at about = 1.963, that is just after 
Q c = 1.9629, where the jump of the angular momentum from L = 8 to L = 12 takes place. 
In Fig. [3] (b) we show a zoom of the relevant part of the figure, where we confirm that this 
structure corresponds to the expected plateau. The total jump in the resistivity is for xq = 0.7 
approximately from p yx = 270 to 291, which is about a 7% of the total value of the resistivity. 
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Figure 4. Hall resistivity p yx as a function of f l when restricting both the current and the 
position operators (see text). Same parameters as in Fig. [3] The presence of a plateau is very 
evident when using the restricted operators. The plateau begins at about f l = 1.963 just after 
= 1.962 where the jump of the angular momentum from L = 8 to L = 12 approximately 
takes place (marked with a vertical gray dotted line). 


We also see that this plateau is soften when the position of the impurity is changed. 

4. Emphasizing the plateau structure 

Up to this point, it has been proved that the most successful way to generate a plateau is by the 
localization of a very special interval of Q where the localized density decreases (see Fig. |2]). 
Within this interval, there is a simultaneous change of B* and p d . The necessary condition 
to produce this effect is under the requirement that for each B* , its ’’dancing pair”, p d , is 
well defined. To be more precise: For each B*, the deduction of p yx requires several steps, 
some of them are operations like ff qp < q\j y \p > or ff qp < q\x\p > where | p > and | q > are 
Fock states: \rii, n 2 , ... ni m ), n* being the orbital occupations. In order to have a fixed value of 
p d the operators j y and x must not do significant changes on the localized part. 

To be sure that this phenomenology is well captured by our model, we forced the 
mechanism in an alternative calculation defining a kind of restricted operators j y and x 
that fulfill the previous requirement: They do not change the localized part of the density. 
We consider in the summations only those elements that couple the vectors with the same 
occupation of the localized orbital (which is n 2 for N = 3 or ri\ for N = 4). The result shown 
in Fig. |4] confirms our intuition. Fig.[4]must be compared with Fig. [3ja). The improvement of 
the visualization of the plateau for some values of xq is evident. 

Within the intervals of linear dependence of p yx with B* where quantum and classical 
behavior coincide llT6l the use of the restricted operators is irrelevant as the occupation of the 
localized orbital is flat. 
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5. Summary and discussion 

Within the framework of quantum simulation, we considered a well known effect for electrons 
under strong magnetic fields. A great amount of experimental data has been accumulated 
in solid state devices (Ref. Ifl6l0 and the most appealing feature lies in the universality of 
the results. The plateaux of the Hall resistivity as a function of the (strong) magnetic field, 
obtained for very clean samples yet containing some disorder, signal the presence of peculiar 
states with well defined fractional filling factors. The values of p yx on the plateaux depend 
only on universal constants (h and e). Up to now, most of these special states still require 
a challenging explanation not always complete. The translation of this physics to finite 
systems of atoms and its interpretation is not always easy. As a first attempt, our goal is 
to understand the mechanism in small samples and try to extrapolate the results to larger 
systems. Notwithstanding the fact that we can only tackle with small number of atoms, our 
results can be experimentally tested, given the ability of the new technologies to deal with 
very small samples and to engineer tunable impurities. 

Summarizing, we have obtained a plateau in the Hall resistivity p yx /Q following the 
expected line of search suggested by the known mechanism in the case of electrons and real 
fields in the fractional quantum Hall regime. The appearance of a flat region on p yx reveals the 
essence of our goal, namely that there is a transfer of atoms from localized to non-localized 
atoms as B* increases. We proved that it is necessary to have part of the system localized 
around an impurity. The transfer from localized to extended states allows the simultaneous 
variation of B* and p d producing a constant value of p yx ~ B* / p d along certain interval of 
Q. Particle interaction and the presence of impurities are crucial ingredients, at odds with the 
case of the integer quantum Hall effect. The improvement of the visualization of the plateau 
is achieved by forcing the complete exclusion of the localized part, using a restricted version 
of the operators for the current and the perturbation. 

So far we have used the expressions involved in the transport equation, valid for 
macroscopic systems with uniform density. Namely, the resistivity (or the conductivity) have 
no space dependence, and the known experimental behavior (see Ref. If54l ) is well captured 
by the linearity of p yx against B* (aside plateau structures or phase transitions of the GS). In 
contrast, our system is finite and the mean density decreases with B* producing a non-linear 
behavior of p yx with B*, which is a finite size effect. 

For systems with larger number of particles (several hundreds or thousands of atoms), 
the presence of impurities represents a small perturbation in the Hamiltonian. Therefore, for 
large systems with impurities randomly distributed, we do not expect that the plateau depends 
on the position of the impurities. This is not the case for systems with a small number of 
atoms, as in this case the presence of the impurities cannot be considered a small perturbation 
and the plateau depends on the position of the impurity (see Fig. [3]). Indeed, the ground state 
is strongly modified and so is p yx , as shown in Fig. 0 with the bump before the plateau in 
Fig.0being due to the presence of the phase transition taking place as the angular momentum 
jumps from L = N(N — 2) to L = N(N — 1). The change of angular momentum modifies 
the resistivity and produces a bump partially overlapped with the plateau. 


Quantum simulation of conductivity plateaux and fractional quantum Hall effect... 

6. Acknowledgments 


12 


We appreciate interesting discussions with Jean Dalibard, Carlos Tejedor, Manel Bosch 
and Bruno Julia-Dfaz. A.T. acknowledges support from the Italian PRIN “Fenomeni 
quantistici collettivi: dai sistemi fortemente correlati ai simulatori quantistici” (PRIN 
2010_2010LLKJBX) and the Universitat de Barcelona for the kind hospitality. We 
acknowledge partial financial support from the DGI (Spain) Grant No. FIS2011-24154 and 
FIS2013-41757-P; the Generalitat de Catalunya Grant No. 2009SGR-1003, 2014SGR940 
and 2009SGR21. We acknowledge also support from EU grants OSYRIS (ERC-2013-AdG 
Grant No. 339106), SIQS (FP7-ICT-2011-9 No. 600645), EU STREP QUIC (H2020- 
FETPROACT- 2014 No. 641122), EQuaM (FP7/2007-2013 Grant No. 323714), Spanish 
Ministry grant FOQUS (FIS2013-46768-P), the Generalitat de Catalunya project 2014 SGR 
874, and Fundacio Cellex. J.T. is supported by grants FPA2013-46570, 2014-SGR-104 and 
Consolider grant CSD2007-00042 (CPAN). Financial support from the Spanish Ministry 
of Economy and Competitiveness, through the Severo Ochoa Programme for Centres of 
Excellence in R&D (SEV-2015-0522) is acknowledged. 


Appendix A. Linear response theory 


Following the standard protocol to analyze the linear response of a system [55], given 
an operator —A Gf(t) simulating a small perturbation (being G time-independent) and an 
operator F for a measurable quantity, the dynamical evolution of the expected value of the 
observable F is given by 

(F)t ~ (F) o = A/(t)|xM|, (A.l) 

where x(w) has been defined as 

<0|F|z/>WG|0) 


xM = El 


E v — Eq + u) + irj 

= |xM|e“<">. 


i /^0 

<0|GHM-F|0} 


(A.2) 


E u — E 0 — lo — irj 

The sum is extended to all excitations and // is a small quantity. 

To obtain information about the Hall response and simulate Eq. ©, connecting with 
Section II, we make the following choice 


N 


G = ^2xi, 

i= 1 

F = Jy, 


(A.3) 
(A. 4) 


and 


f{t) = £(t)sin(u}t), 


(A.5) 
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Figure Al. (a) Comparison between the calculation of the Hall resistivity p yx as a function 
of fl with the complete dynamical evolution (TE) and the result from the LRT for N = 3 and 
ai = (1, 0). We use here the restricted operators. Notice that for N = 3 within this range 
of O no plateaux occur, (b) same fr N = 4 when a plateau is apparent. The figure shows that 
even in this case, we obtain nearly the same result. 



Figure A2. Hall resistivity p yx as a function of SI for N = 3 and ai = (1,0) . This figure 
shows a typical resonant structure in the LRT calculation (upper curve) which disappears in 
the TE case (lower curve). 


see Eqs. (fl4l) . (fl5l) . and (fT8l) . If we identify the perturbation with the electric potential 
associated with a space independent electric field directed along the X-axis, 

E x (t) = Xf(t) (A.6) 

then, from Eq. (IA.1I) it is 

<Tyx = IxMl ( A - 7 ) 

and from Eq. (l20l) we obtain the transverse resistivity p yx . The function x( u ) obtained 
numerically fulfills the condition x*(x) — x (— u ) as wc h as the Kramers-Kronig 

relations lf55l . 

In what follows we consider three significant results. Fig. lAll faf shows a comparison 
between the two techniques: TE and the LRT for p yx /Q. This result indicates that, for a wide 
range of the Hall effect is a linear response process. Figure I A1 1 (b) is again a verification 
that the Hall response is linear even in the case where the response implies the formation of a 
plateau. Finally, we include Fig. IA2l to illustrate possible limitations of the LRT. Notice that 
the structure of the upper part of Fig. IA2I is the typical result of a resonance: For a critical 
value of and for some excitations, a denominator in the expression of ;\;(co') [see Eq. (I A. 21) 1 
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nearly vanishes. The system absorbs energy and the behavior is not linear. Namely the LRT is 
not applicable there. It could well be a sign of a peculiar state and its Hall response. However, 
the lower part of the figure demonstrates that this is not the case, that it is not related with the 
phenomena we are looking for. Notice that Fig. lA 1 1 (a) is for N — 3, in this case, p yx /B* does 
not present a plateau, at odds with the result shown in Fig. lAll (b) for N = 4 where within the 
same range of Q values, a plateau is apparent. In Fig. IA11 the LRT result does not show the 
resonance due to the low resolution. 
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